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Abstract. We study dynamical aspects of circumstellar environment around massive zero-metallicity first stars. For this pur- 
pose we apply our NLTE wind models. We show that the hydrogen-helium stellar wind from stationary massive first generation 
(Population III) stars (driven either by the line (bound-bound) or continuum (bound-free and free-free) transitions) is unlikely. 
The possibility of expulsion of chemically homogeneous wind and the role of minor isotopes are also discussed. Finally, we 
estimate the importance of hydrogen and helium lines for shutting off the initial accretion onto first stars and its influence on 
initial mass function of first stars. 
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1. Introduction 

The Big Bang nucleosynthesis has left our universe completely 
metal-free (e.g. Coc et al. 120041 . Thus, at the dawn of the 
universe there existed a population of stars composed from 
hydrogen and helium only. Numerical simulations show that 
the physical properties and the evolution of these first stars 
(Population III stars) were very different from the stellar evo- 
lution and properties of Population I and II stars. Simulations 
of first star formation (cf. Bromm et al. [1999 Nakamura & 
Umemura 120021 show that due to the less efficient cooling 
the first stars were probably much more massive than the 
present stars. First stars with masses of the order 100 Mq might 
have been formed. However, there are indications that the ini- 
tial mass function of first stars was bimodal (Nakamura & 
Umemura 2002 Omukai & Yoshii 2003 1, and thus solar mass 
first stars may have been formed as well. 

There are no available observations of first stars up to now. 
Although the recently discovered extremely low metallicity star 
HE 0107-5240 (Christlieb et al. 120021 was claimed to be the 
first star (Shigeyama et al. 2003 1, there are some problems with 
this explanation. It is difficult to explain the present surface 
chemical composition of this star (especially that of carbon and 
nitrogen) on the basis of initial pure hydrogen and helium com- 
position. For example, Shigeyama et al. (2003 1 explain the CN 
composition by the surface enrichment due to the central nu- 
cleosynthesis, however there is a question whether significant 
amount of CN elements may occur on the surface for this stel- 
lar type according to evolutionary calculations of Marigo et al. 



2001 1. It is more likely that this star was formed from the ma- 
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terial which was already enriched by the first stars (Bromm & 
Larson 2004 1 or by subsequent stellar populations. 

Consequently, it is necessary to constraint ourselves only 
to indirect observation of the first stars. Probably the most re- 
liable indirect observation of the first stars comes from the po- 
larisation measurement of cosmic microwave background radi- 
ation. If the first stars ionized the matter in the universe sig- 
nificantly, we should observe its imprint on the cosmic mi- 
crowave background due to the light scattering on free elec- 
trons. Recent cosmic microwave background polarisation mea- 
surement showed that the universe was reionized 180 x 10 6 
years after the Big Bang (Kogut et al. 2003 1. Detailed numeri- 
cal analysis shows that this effect can be attributed to the gen- 
eration of massive zero-metallicity stars (Cen 2003 Wyithe & 
Loeb 2003 Sokasian et al. 2003 1. Moreover, the level of the 
sky infrared background radiation can also be explained by the 
generation of massive first-generation stars (Magliocchetti et 
al. 120031 . 

The detailed theoretical study of physical properties of the 
first generation stars is a challenging task. Especially appealing 
is the study of its stellar winds, since winds might have a sig- 
nificant impact on the first stars evolution. Moreover, the study 
of outflows of the first stars may be important for the estimate 
of the duration of supermassive stars formation. The possibil- 
ity of supermassive stars formation depends critically on the 
metal abundance in the interstellar medium. The presence of a 
relative small amount of metals might inhibit subsequent for- 
mation of supermassive stars (see Bromm & Larson 2004] for 
a review). Note, for example, that there are indications that the 
initial formation of supermassive stars might have been inhib- 
ited by metals expelled due to the supernova explosions and 
that direct black hole formation by implosion is necessary for 
the explanation of the duration of massive star era (Ricotti & 
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Ostriker 2003 1. Thus, the interstellar medium had to remain 
metal free for a significant amount of time if the stellar wind 
did not enrich it with nuclearly processed material. This shows 
how important the study of winds from first generation stars 
potentially could be. 

The outflows from first stars have not been systematically 
studied up to now. First step in this direction was done by 



Kudritzki ( 2002 1, who calculated wind parameters suitable for 
massive extremely-low metallicity stars. As was discussed by 
Kudritzki and also subsequently supported by detailed numer- 
ical modelling by Krticka et al. (|2003 ), some of these stellar 
winds may suffer from multicomponent effects (see Krticka & 
Kubat (2001 1 for a detailed modelling of multicomponent ef- 
fects in the line-driven stellar winds). However, to our knowl- 
edge the line-driven stellar winds of zero-metallicity stars have 
not been discussed so far. To fill this gap we present the detailed 
discussion of pure H-Hc line-driven winds. 

2. Theoretical limits for line-driven winds 

The necessary condition for launching the line-driven stellar 
wind is relatively simple. The radiative force should exceed the 
gravity force at some point. Thus, instead of presenting more 
detailed numerical models we start with a discussion of this 
condition and its application to simple situations. This helps 
to better understand more complicated behaviour of numerical 
models. 



2. 1. Minimum radiative force 

The radiative acceleration due to the absorption or scattering of 
radiation is given by (Mihalas lll978> 

A /-CO 

ff rad = — / XuH v dp, (1) 



cp Jo 

where p is the mass density, \ v i s me absorption coefficient, 
H v is the radiative flux, and v is the frequency. For the calcula- 
tion of the radiative force in the optically thin environment, the 
radiative flux H c (v) emerging from the stellar atmosphere can 
be directly inserted into Eq. Q. However, this is not the case 
of the stellar wind since the radiative flux H v is influenced by 
the absorption and emission of the radiation in the sorround- 
ing wind environment in such case. Thus, the self-shadowing 
in the given transition has to be accounted for, preferably by 
the solution of the radiative transfer equation. In the case of 
rapidly accelerating stellar wind a convenient approximation of 
the radiative transfer equation solution is given by the Sobolev 
approximation (see Sobolev [l947> . 

From Eq. Q follows the limit below which classical line- 
driven wind is not possible. In particular, neglecting continuum 
absorption, for a given set of lines the radiative force is max- 
imum, if lines are not self-shadowed and if lines are optically 
thin, i.e. if the optical depths t < 1 in lines (Gayley Q995J. 
Thus, for given occupation numbers the line acceleration (i.e. 
line force per unit of mass) cannot be larger than 

47r 2 e 2 
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where now H c (vij) is the frequency dependent flux emerging 
from the star at the stellar radius i?,*, (l — /j 2 ) = R 2 /r 2 is a 
correction due to the dilution of radiation, z/y and gifij are the 
frequency and the oscillator strength of a given transition i *-* 
j, m, rij are the number densities, and gi, gj are the statistical 
weights of corresponding levels. 

Note that g rad = max i s related to the Gayley's parameter Q 
(Gavley [T995l Puis et al. l2000l namely 
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where <?£ ad is the radiative acceleration on free electrons, 
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where L is the stellar luminosity, 
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s e is the Thomson scattering cross-section, and n c is the elec- 
tron density. The line-driven wind is possible only if the total 
radiative force is greater than gravity at some point in the stellar 
atmosphere or circumstellar envelope. Thus, the necessary con- 
dition for line-driven wind reads (neglecting the gas-pressure) 

^rad, max + ff rad = (Q + l) > ( 6 ) 

where g = GM/r 2 is the gravity acceleration. Using Eq. 10 
for maximal line acceleration and 1 — p 2 = R 2 jr 2 we can 
obtain a condition for the lowest limit of line-driven wind as 
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where M is stellar mass and 
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Using the Gayley's parameter, which for our case reads 
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where 7J t (i?„) = i/(47ri?») 2 is the total flux at the stellar 
surface i?*, the condition l|6} may be rewritten as 
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(10) 
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Let us first assume that the stellar surface flux H c {vij) is con- 
stant for a given transition i — » j throughout the wind, i.e. that 
the frequency dependence of H c (vij) due to the Doppler shift 
can be neglected. In such a case the only quantities in Eq. 
which depend on radius are the fractions fii/p, nj/p. Thus, for 
a given star Eq. Q gives a lower limit for rij j p below which 
the line-driven wind is not possible. Eq. is a necessary con- 
dition for launching the stellar wind. 

Another formulation of the condition for launching the stel- 
lar wind may be obtained from the requirement that the work 
done by the radiative force (per unit mass) is greater than the 
absolute value of gravitational potential. It is noteworthy that in 
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the case of wind driven purely by optically thin lines this con- 
dition can be easily obtained by the integration of Eq. (|6j from 
r = i?* to infinity if the wind ionization and excitation state 
is constant (note that in this case fraction rij/p is constant for 
homogeneous winds) 



YGM 
I! 
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(ii) 

Apparently, this is an identical condition to Eq. and thus 
Eq. is also an energy condition for launching the stellar 
wind driven purely by optically thin lines. However, real line- 
driven stellar winds are accelerated not only by optically thin 
lines but by the mixture of optically thin and optically thick 
lines. Thus, a simple condition (II It cannot be used in this case. 
Moreover, some of the lines which accelerate the stellar wind 
may be optically thick in the photosphere. In such a case, the 
stellar flux H c depends on the wind velocity (and thus also on 
the radius) due to the Doppler effect and condition ( II It cannot 
be applied. 

2.2. Application to hydrogen 

The necessary condition for the existence of line-driven winds 
is especially simple in the case of the radiative force only 
due to the hydrogen lines. This is only a rough approximation 
even in the case of stars with zero metallicity for which also 
the radiative force due to the helium lines may be important. 
However, condition at least provides a guess of occupation 
numbers for which the line-driven winds of stars at zero metal- 
licity may exist. Moreover, adding other hydrogen-like ions is 
straightforward. 

The oscillator strengths for hydrogen are given by the well 
known Kramers formula, which neglecting Gaunt factors reads 
(Mihalas fT978l Eq. (4.78) therein) 



/V mm — 
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where mn is hydrogen atom mass, X — n-nm-n/ p is the 
relative hydrogen abundance, nn is hydrogen number density 
and N^l\ = n^f\/nji (where n^"\ is the minimum number 
density of ground state neutral hydrogen necessary to launch 
the stellar wind). Apparently, for stars close to the Eddington 
limit (T — > 1) the minimum relative number density Nm.i for 
launching the hydrogen line-driven wind is lower. 

2.3. Application to helium 

He II is the most important helium ion for radiative driving 
since helium is ionized in the stellar winds of hot stars. Thus, 
the radiative force due to helium can be calculated in a very 
similar way to that due to H I. Similarly to the case of H I, the 
resonance lines are the only important He II lines for radiative 
driving (assuming that helium is singly ionized ionized). Thus, 
we can write for He II (similarly as for H I) 



9i 
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The condition of existence of wind driven by the absorption in 
He II lines is (see Ea.ll4> 
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and the lowest relative occupation number iVHen,i f° r which 
the helium-driven wind is possible is 
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For the resonance lines it may be assumed 
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(12) 
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which in fact means neglecting the stimulated emission in all 
resonance transitions. Moreover, the contribution of other than 
resonance lines to the radiative force can be in many cases ne- 
glected. 

Thus, in the case of hydrogen, we can obtain from Eq. 
relatively simple condition for the existence of line-driven wind 
in the form of 
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(14) 

The main contribution in the sum in condition J14t is due to 
the resonance lines and the lowest relative occupation number 
of the hydrogen ground level A^hi.i, for which the hydrogen- 
driven wind is possible, is given by 



iv HeIl,l 



-3 



3V3m c m H cC 2 GM 
128ne 2 YR 2 



(1-T), (18) 



where mnc is helium atom mass, Y = tt-hc?71Hc//0 i s tne re l" 
ative helium abundance, hhc is helium number density and 
-^Hen l = n Heu i/ n He- The sums in H7\ and J 1 81 converge 
due to the asymptotic behaviour of hydrogen-like oscillator 
strengths, gifij ~ 1/j 3 (see Eq. H2\ ). 

In many cases the stellar wind is accelerated by lines which 
are optically thick in the photosphere and thus the stellar flux 
H c for frequencies close to the photospheric line center de- 
pends strongly on frequency. Since the radiative flux at a given 
point is shifted due to the Doppler effect, the optically thin ra- 
diative force depends on velocity (and consequently also on 
radius). Even in such a case the conditions ( I15I18> can also be 
used to estimate the minimum conditions for line driven winds. 
However, the actual radiative flux H c [v(v r )], which depends 
on wind velocity (and thus also on radius) should be inserted in 
these equations. 
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2.4. Stellar wind driven by optically thin lines 

We have derived the necessary condition for launching the line- 
driven stellar wind. However, there is another question which 
should be answered. Is this condition also the sufficient one? 
In other words, are the optically thin lines only able to drive 
stellar wind? 

2.4.1. Constant ionization fraction 

Let us first assume that the wind ionization and excitation state 
does not depend on wind density and on radius and that the stel- 
lar surface flux H c does not significantly vary with frequency. 
In such a case the optically thin radiative acceleration Q does 
not depend on wind density. Thus, also the wind momentum 
equation 

dv r 2a 2 a 2 dv r thin d 

v r -7— = 1 j 9 + 9 +9 C , (19) 

dr r v r dr 

(where the optically thin radiative acceleration g thln = 
<7 rad ' max , a is sound speed, and where we assumed isothermal 
flow) does not depend on wind density. It may seem that the 
optically thin line-force may launch the stellar wind with an ar- 
bitrary mass-loss rate. However, this is not true, since in such 
a case either for higher wind densities some of the lines would 
become optically thick or the photon tiring effect (Owocki & 
Gavlev lT997> would impose constrains on the mass-loss rate. 

2.4.2. Variable ionization fraction 

In the case when the wind ionization and excitation state de- 
pends on wind density, the value of wind mass-loss rate can be 
derived from the hydrodynamic equations. The optically thin 
momentum equation dl9l has a critical point for v r = a. To 
obtain a smooth wind solution, the regularity condition 

2n 2 

— +g thin + g? d = 9 (20) 

r 

must be fulfilled at the critical point. Since wind ionization and 
excitation state depends on wind density now, also the optically 
thin radiative acceleration depends on wind density. Thus, the 
critical condition i20\ depends implicitly on the wind density 
and using this equation the mass-loss rate of thin winds can be 
calculated. Note that in the common case when the gas pressure 
term is negligible, the condition will be just fulfilled at the 
critical point. This again shows the importance of the condition 

However, this analysis does not tell us anything about the 
possibility to drive the wind material out from the stellar grav- 
itational potential well. This condition has to be tested using 
detailed numerical modelling. On the other hand, more wind 
properties can be unveiled using a simple analysis. The stellar 
photosphere deep below the wind region shall be quasi-static. 
Thus, to obtain a smooth flow, the radiative force has to be 
smaller than the gravity in this region and the condition has 
not to be fulfilled below the critical point. The change of wind 
ionization or excitation state then has to assure that the equality 
in (|6j is fulfilled at the critical point. Finally, in the outer parts 



of the wind the inequality in condition has to hold, i.e. the 
radiative force has to be greater than the gravity in the outer 
wind parts. 

2.4.3. Frequency dependent stellar flux 

The frequency variable stellar flux gives probably the most at- 
tractive possibility to obtain consistent radiative force to drive 
an optically thin wind, i.e. radiative force which is lower than 
the gravity force below the critical point and higher outwards. 
This is caused by the fact that some lines, which accelerate the 
stellar wind, may be optically thick in the photosphere and op- 
tically thin in the stellar wind. In such a case the stellar flux 
H c strongly depends on frequency and, consequently, due to 
the Doppler effect, also on the wind velocity. The stellar wind 
near the star has low velocity and a given line is accelerated 
mainly by the weak flux from the line core. On the other hand, 
stellar wind with higher velocities is accelerated mainly by the 
stronger flux from the line wings and thus the line force of the 
outer wind may be higher than that of inner part of the wind 
close to the star. However, this is not as simple as it looks like 
since higher flux may also cause higher excitation and ioniza- 
tion and consequently lower the radiative force. Clearly, de- 
tailed numerical calculations are necessary to test this possibil- 
ity. 

3. Numerical calculation of thin winds 

For the numerical tests we used NLTE wind models described 
by Krticka & Kubat (2004i. The computer program allows to 
solve hydrodynamic, simplified radiative transfer, and statis- 
tical equilibrium equations in a radiatively driven stellar wind. 
However, for the present purpose we slightly changed our code. 
Since we do not know in advance whether the stellar wind is 
possible or not, the hydrodynamical variables, i.e. velocities, 
temperatures, and densities of all wind components (with an 
exception of the electron density) are kept fixed. This enables 
us to calculate the model occupation numbers and the radiative 
force regardless the existence of the wind and, consequently, to 
test the wind existence. 

The method itself was slightly improved with respect to 
the published one. We included accelerated lambda iterations 
(cf. Rybicki & Hummer 1992 1 and Ng acceleration (Ng |1974> . 
These improvements will be described in detail in a separate 
paper (Krticka & Kubat, in preparation). Moreover, to obtain a 
correct surface flux we accounted for the Doppler effect dur- 
ing the calculation of radiative transfer in lines, i.e. we shifted 
the stellar surface flux H c according to the actual wind veloc- 
ity (we accounted for the detailed profiles of photospheric lines, 
cf. Babel 1996). Finally, to assure a smooth transitions between 
stellar atmosphere and wind, the model atmosphere fluxes are 
taken at the point where the density is equal to the boundary 
wind density. 

First, we have performed tests whether the radiative force 
is large enough to drive stellar wind in optically thin lines. 
Although this gives only the necessary condition for stellar 
wind and not the sufficient one (and thus further calculations 
are necessary to test whether the stellar wind is possible or 
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not), this approach enables us to better understand more de- 
tailed models. 

3. 1. Minimum number density for line-driven winds 

We performed two different tests of the possibility of a stel- 
lar wind of zero-metallicity stars driven by optically thin lines. 
First, we calculated the ionization structure of a possible ex- 
tended stellar atmosphere and compared this ionization struc- 
ture with minimum number densities for launching a line- 
driven wind obtained using Eqs. (1151 and Jl 8I >. 



Table 1. Stellar and wind parameters used to calculate dis- 
cussed wind models taken from Kudritzki (2002). 





Stellar parameters 




Model 


log (L/Lq) 


Tcff 


log 5 


R, 


Twind 






[K] 


CGS 


[R©] 


[K] 




7.03 


60 000 


3.95 


30.38 


49 000 


T60 




50 000 


3.63 


43.76 


39 000 


T50 




40 000 


3.25 


68.32 


29 000 


T40 


6.91 


60 000 


3.99 


26.24 


49 000 






50 000 


3.68 


38.06 


39 000 






40 000 


3.28 


59.49 


29 000 




6.76 


60 000 


4.04 


22.29 


49 000 






50 000 


3.73 


32.10 


39 000 






40 000 


3.34 


50.12 


29 000 




6.57 


60 000 


4.11 


17.78 


49 000 






50 000 


3.79 


25.74 


39 000 






40 000 


3.41 


40.18 


29 000 




6.42 


60 000 


4.16 


15.07 


49 000 






50 000 


3.85 


21.71 


39 000 






40 000 


3.46 


33.93 


29 000 




6.30 


60 000 


4.21 


13.11 


49 000 






50 000 


3.89 


18.88 


39 000 






40 000 


3.50 


29.51 


29 000 





For our detailed discussion we selected only those stars 
from the Kudritzki (2002i list, which have the largest lumi- 
nosity, i.e. those which are close to the Eddington limit. These 
stars have the largest radiative force due to the light scattering 
on free electrons in their atmospheres and thus also the highest 
chance to launch the stellar wind. Detailed models were calcu- 
lated also for other stars from the Kudritzki's list, however for 
these stars we describe the final results only. Parameters of all 
studied model stars are given in Tab.[2 

We compared the minimum number densities necessary for 
launching the stellar wind with actual number densities calcu- 
lated by our NLTE code. Since we do not know in advance 
whether the radiative force is strong enough to drive a wind, 
the velocity structure of our models is given by an artificial ve- 
locity law 

?; ( r ) = 10-3 /^£ +4 10 7 cms- 1 ^— (21) 
V 3 m-H R* 

This velocity law corresponds to the Taylor expansion of the 
wind velocity around the sonic point. According to our expe- 
rience with wind models of O stars this expression is valid in 



the wind region close to the star with r < 1.2i?„ within the 
accuracy of 20%. This is the region of rapid wind accelera- 
tion where the stellar wind attains velocities of the order of 
100 km s -1 . Note, however, that obtained results do not signif- 
icantly depend on the detailed form of the velocity law. The 
density structure is obtained from the equation of continuity. 
Wind temperature in these models was equal to T w j n d and elec- 
tron density was consistently calculated from hydrogen and he- 
lium ionization balance. For each of these stars with differ- 
ent effective temperatures three models with different values 
of mass-loss rates were calculated, namely 10~ 6 Mq year -1 , 
1O -7 M year -1 , 1O -8 M year -1 . 

Comparison of model number densities and minimum 
number densities necessary to launch the wind is given in 
Figs. [0-0 Near the stellar surface number densities of ground 
levels of hydrogen-like ions are relatively high and they de- 
crease with radius due to increasing ionization. The ionization 
is higher in the outer regions due to both lower density and 
higher radiative flux at frequencies outside the line core (due to 
the Doppler effect the line frequency shifts from the saturated 
core of photospheric line-profile, consequently larger flux ab- 
sorbed in the lines increases populations of higher levels, for 
which the ionization rate is much higher than for lower levels). 
Minimum number densities also decrease. This is caused by in- 
creasing stellar flux as the wind lines move out from the centre 
of photospheric line due to the Doppler effect. 

The relative number density of the ground level of H I is 
very low, mostly much lower than the minimum number den- 
sity necessary to launch the stellar wind by H I lines. Thus, we 
conclude that H I lines are not important for line-driven winds 
of first stars. On the other hand, for specific regions of some 
models the actual number density of ground level of He II is 
higher than the minimum number density necessary to launch 
the stellar wind. In these regions the total radiative force may 
be higher than the gravity force and the sufficient condition for 
launching the stellar wind may be fulfilled there. 

3.2. Maximum radiative force 

The preceding simple analysis can be extended using the de- 
tailed calculation of the radiative force. For stars in Tab.^we 
have calculated the relative difference (,g thm + g l c ad — g) /g 
(i.e. the net acceleration in terms of gravity acceleration) as 
a function of radius. Corresponding graphs are plotted in 
Figs.S-i] 

We came to the same conclusion as in the previous para- 
graph. The maximum possible radiative force is greater than 
the gravity in some regions of some models. However, the func- 
tional dependence of the difference g thm + g£ ad — g in most 
cases does not have a correct radial dependence necessary to 
launch the stellar stellar wind, i.e. the total radiative force in 
these models is higher than the gravity near the star and lower 
than the gravity in the outer regions. 

Moreover, we have to keep in mind that the optically 
thin radiative force is maximum for a given set of lines and 
clearly the actual radiative force will be lower due to the self- 
shadowing effect. Moreover, for many of the models the total 
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Fig. 1. Comparison of minimal occupation numbers -/V^ 1 ", 
N^l™ 1 required to launch the wind and actual occupation 
numbers iVm i and AT HeIIjl (shown by unlabeled curves) for 
model T40 and different mass-loss rates of 10~ 6 Mq year -1 , 
1CT 7 M Q year -1 , and 10" 8 M Q year -1 . 
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Fig. 2. The same as Fig. [2 however for model T50. 

Inspecting Figs. [Q-[6]note that Eqs. (I15> and Jl 8i are able 
to reliably estimate the minimum occupation number neces- 
sary to launch a line-driven stellar wind because they correctly 
predict the location of region where the optically thin radiative 
force is greater than gravity. 



radiative force is not greater than the gravity. Thus, we con- 
clude that winds of zero metallicity stars driven purely by the 
optically thin lines are unlikely. However, these results shall be 
later firmed up by a more detailed models with the possibility 
that some lines are optically thick. 



4. Wind with optically thick lines 

Our final test of the possibility of the existence of line-driven 
stellar wind of first stars should include the natural possibility 
that some lines are optically thick. Thus, we now turn our atten- 
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Fig. 3. The same as Fig.^ however for model T60. 



tion to more realistic case of the stellar wind driven partially by 
the optically thick lines. This is a much more complicated task 
than the wind driven purely by optically thin lines since now we 
cannot use any simple criterion to decide whether the wind is 
possible or not. We have to use detailed modelling. Moreover, 
the nonlinear character of wind equations (especially equations 
of statistical equilibrium) does not guarantee that the obtained 
solutions are unique. This is especially the case of zero metal- 
licity stars when the wind may be driven just by one ionization 



state. Simply, for some ionization structures the wind may be 
possible and for some other not. 

To make this problem more tractable we calculate radia- 
tive force in the Sobolev approximation for the same models, 
for which we studied optically thin radiative force in Sect.|3| 
The radiative force in the Sobolev approximation is given by 
the sum of partial radiative forces due to the individual lines 
(Castor [T974l Abbott Il982l 



^rad 



■iTT V r ST^ tt 

9 =—^—y.VijH c 

r lines J ^ c 



/i (l + crpL 2 ) (l - e T ») d/i, 



(22) 

where p is the wind density, v r is the radial velocity, H c is the 
stellar flux, \i c = -Jl — R 2 /r 2 is the direction cosine of the 
photospheric limb, i?» is stellar radius, a = dlnw r /dln?' — 1 
and the Sobolev optical depth r M is 



7re 



-) Si/i 

9i 9j 



13 v r (1 + a^) ■ 



(23) 



If some of the lines are optically thick, then the actual radia- 
tive acceleration given in the Sobolev approximation Eq. i22\ 
is lower than the maximum radiative acceleration (|2j. This is 
caused by the effect of self-shadowing, i.e. that some part of 
the stellar flux which is supposed to accelerate the wind in a 
line was already absorbed by a given line in the region closer 
to the star. Thus, here we present only results for models for 
which optically thin radiative force can drive a stellar wind, 
since it is clear that there will not be any wind in the opposite 
case. 

The radiative force calculated using Eq. &22\ in the selected 
models is given in Fig.0 Clearly, since most of the strongest 
lines are optically thick, the radiative force due to these lines 
is much lower than the optically thin radiative force. More re- 
alistic radiative force given by Eq. (1221 is generally by several 
orders of magnitude lower than the gravity force. Thus it is un- 
likely that some change of wind density or velocity may give 
significantly stronger radiative force which would drive a stel- 
lar wind. This is also supported by our models, because we 
were not able to find sufficiently large radiative force for any 
reasonable wind models. We conclude that it is unlikely that 
stationary massive H-He stars have any H-He stellar wind 
driven by the line transitions. 

We have also calculated similar models for other stars from 
Kudritzki ( 2002 ) list with lower masses and lower luminosities 
(with parameters given in Tab.Q. These stars are not so close to 
the Eddington limit as stars discussed sofar. Consequently, the 
radiative force due to light scattering on free electrons is lower 
for these additional stars and the minimum radiative force nec- 
essary to launch the stellar wind is higher in this case (see 
Eas. ll5l[T8l . Clearly, the possibility for launching of the stellar 
wind is lower for stars which are not so close to the Eddington 
limit. 

5. Influence of continuum radiative force 

Another force which can potentially accelerate the stellar wind 
of first stars is the radiative force due to bound-free and free- 
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Fig. 7. The plot of the net radiative acceleration g 1 ' + g' c — g and the net radiative acceleration of the whole gas due to free 
electrons only gj ad — g relative to the gravitational acceleration g. Clearly, the contribution of the radiative acceleration due to 
the line transitions g rad calculated after Eq. (I22t is much smaller than the gravity acceleration and is not able to drive a stellar 
wind. 



free transitions. Numerical calculation showed that the contri- 
bution of free-free transitions to the radiative force in the cir- 
cumstellar environment is much lower than that of bound-free 
transitions. Thus, we shall discuss the bound-free transitions 
only, although the free-free transitions are also included in our 
models. The radiative acceleration has in this case the form 
(Mihalas [T978l 



„rad,bf . 



47T 



Xv 



(24) 



where \v is the opacity due to the bound-free transitions, 



a ik {v), (25) 



where summation goes over all levels i for which bound-free 
transitions are accounted for and the asterisk refers to the LTE 
value. 

Let us first briefly discuss an interesting possibility that 
the stellar wind is driven purely by the bound-free absorp- 
tion. Combining Eqs. d24l and d25l . the radiative acceleration 



is given by 
9^ = - 



P 



a lk (v)H v e, 



dv 



(26) 



As has been already mentioned in Sect. 12.41 in order to obtain 
smooth wind solution, the radiative force has to be smaller than 
the gravity force in the stellar atmosphere below the critical 
point and larger than the gravity force above the critical point. 
It is very difficult to achieve this in this case, since the radiative 
flux H v decreases with radius rapidly, especially in the case 
when strong bound-free absorption is present. Moreover, the 
critical condition is not a sufficient one to drive the stellar wind, 
since it does not tell us anything whether the radiative force in 
the outer wind regions is strong enough to lift the wind material 
out from the stellar gravitational potential well. The decrease of 
the stellar flux with radius due to the bound-free absorption (see 
also the photon tiring effect, Owocki & Gayley 1997 1 can com- 
plicate this. Note that the situation with the line-driven wind is 
simpler, because the lines are usually illuminated by the unat- 
tenuated stellar flux thanks to the Doppler effect. 
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We included the radiative force due to the bound-free and 
free-free transitions Eq. (I24> into our models. The numerical 
calculations (for stars with parameters given in Tab.^ showed 
that the radiative force due to the bound-free and free-free tran- 
sitions does not significantly contribute to the acceleration of 
stellar wind of first stars because it is typically at least two or- 
ders of magnitude lower than the gravity force. 

Finally, we calculated also static spherically symmetric 
NLTE model H-Hc atmospheres (for the description of the 
code see Kubat 2003 1, where the radiative force was calculated 
without any simplification. It turned out that all of the models 
with parameters from the Tableware sub-Eddington i.e. the to- 
tal radiative force is always lower than gravity. However, for 
the most luminous models the value of T reached the value of 

0. 99. 

6. Expulsion of individual elements by a wind 

Whereas first stars probably fail to drive H-He outflow, there 
is a possibility that hydrogen or helium are expelled separately, 

1. e. that there exists a pure hydrogen or helium wind. To show 
whether such outflow of chemically pristine elements is possi- 
ble or not, we recalculated models with the same stellar param- 
eters presented in the Sect.[3] assuming that one element may 
form the wind while the other not and compared the radiative 
force in such chemically homogenous wind with the gravity 
force acting against such wind. To avoid charging of the star 
we assume that a corresponding number of electrons follow 
the escaping element making the wind electrically neutral. We 
assume that such wind exists above the H-He photosphere. 

Since these models are chemically homogenous, the elec- 
tron density in such models is calculated only from the con- 
tribution of the element that remained in the wind (i.e. either 
hydrogen or helium). This difference is crucial for the final 
results, at least in the case of helium. If we assume that the 
helium gas consists of only the isotope 4 He, we may say that 
helium has two times less electrons per nucleon (proton or neu- 
tron) than hydrogen. Hence, the radiative force per unit of mass 
(i.e. the radiative acceleration) due to free electrons is lower in 
the case of helium. Consequently, although the number of free 
electrons per atomic nucleus is higher in pure He than in H-He 
models, the total radiative force is lower. Surprisingly, although 
helium lines significantly contribute to the radiative force in the 
H-Hc models, the possibility of pure He outflow is lower than 
that of H-He. 

On the other hand, also the situation in pure hydrogen mod- 
els is different. Hydrogen plasma has higher number of free 
electrons per nucleon than H-He plasma. Consequently, H-He 
stars which are close to the Eddington limit may have slightly 
higher radiative force than gravity force for the same stellar pa- 
rameters (mass, radius, luminosity) in a pure hydrogen wind, 
whereas the radiative force in hydrogen-helium atmosphere is 
(due to the contribution of heavier helium atom) lower than 
gravity force. 

These considerations are supported by our numerical mod- 
els. The total radiative force in helium models is lower than 
gravitational force in all models considered (for stars with pa- 
rameters given in Tab.[0- On the other hand, for model stars 



which are very close to the Eddington limit, the radiative force 
in purely hydrogen models may be greater than gravity, thus 
potentially enabling hydrogen wind. 

Condition for the existence of such hydrogen wind is rel- 
atively simple. If the wind exists, radiative acceleration due to 
the free electrons in a pure hydrogen plasma is greater than the 
gravitational acceleration, 



7e 3 H > 9, 



(27) 



or using Eddington parameter Th in a pure hydrogen plasma 
(cf.Eq.03 

Hr 

r H = ; > 1, (28) 



where 



4ncGM 



PH 



(29) 



s e is the Thomson scattering cross-section, and is the elec- 
tron number density in hydrogen plasma with the density pu- 
Condition ( 128 1 can be rewritten in a more convenient form. 
Now we compare conditions in pure hydrogen wind with that 
in hydrogen-helium atmosphere. Using Eqs. (TJ} and Q for 
the hydrogen-helium atmosphere, we can express the fraction 
Lj (iwcGM) in fl28i using the Eddington parameter in the stel- 
lar atmosphere T and we obtain 



r = r. 



On.) 



atmosphere 



(^c H )wi 



(30) 



wind 



which yields the condition (p = p H + p He is the total density in 
the atmosphere) 



r>|*l 



wind V ' PHe / atmosphere 



(31) 



For a chemical composition given by primordial nucleosynthe- 
sis and assuming a fully ionized gas both in the atmosphere and 
in the wind we obtain simple condition 



T > 0.859. 



(32) 



Clearly, pure hydrogen-helium stars close to the Eddington 
limit (with 0.859 < T < 1) may have a pure hydrogen wind, 
but not the H-He wind. 

The only stars from our sample for which T > 0.859 and 
which thus can have a pure hydrogen wind are the most massive 
stars with M — 300 Mq (for all considered values of effective 
temperature) and star with M = 250 Mq and T eff = 40 000 K. 

However, despite this behaviour of the radiative force in 
pure hydrogen models, pure hydrogen winds need not necesser- 
ily exist in all stars for which they are theoretically possible. 
Radiative force in these models is not due to the hydrogen line 
transitions, but mostly due to the free electrons. Thus, in or- 
der to hydrogen wind exists, there should be some proccess 
which is able to separate hydrogen and helium atoms. The 
mass-loss rate of such a theoretical wind depends mainly on 
the rate of this separating process. The only possible process, 
which comes into consideration, is the gravitational settling due 
to the different atomic masses (see Michaud 2005 for a recent 
review). However, the typical rate of such a process is very 
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small, thus mass-loss rates of a pure hydrogen wind would be 
probably also very small (see also AppendixfA). Moreover, any 
macroscopic motion in the stellar atmosphere may inhibit this 
gravitational settling. Thus, we conclude that any possible pure 
hydrogen wind has either a very low mass-loss rate or is even 
missing. 

7. Other isotopes and elements 

7. 1. Contribution to the radiative force 

Hydrogen and helium isotopes 2 H and \ He were produced in 
a nonnegligible amount during the era of the Big-Bang nucle- 
osynthesis. The positions of their lines are slightly shifted with 
respect to those of JH and |He and thus they can potentially be 
exposed to a slightly higher flux than their much more abundant 
counterparts. 

To test their importance we included 2 H and jHe into our 
atomic list. We used basically the same atomic data for the cal- 
culation of model atoms, however with modified energy levels 
due to the isotopic shift. The isotopic abundances were taken 
from the simulation of primordial nucleosynthesis by Coc et 
al. ©04 1, 7V(jH)/./V(?H) = 2.6 x 10 -5 , N (JH) / TV (|He) = 
1.04 x 10 -5 , which are in a relatively good agreement with the 
observed values. 

The importance of other hydrogen and helium isotopes was 
tested using optically thin models discussed in Sect.[3] The cal- 
culated excitation and ionization states of 2 H and ^Hc are very 
similar to those of }H and |Hc. Moreover, due to the rela- 
tively small isotopic shifts and due to the relatively broad pho- 
tospheric lines the flux at the positions of lines of less abundant 
isotopes is similar to that of more abundant isotopes. Finally, 
the abundance of these isotopes is several orders of magni- 
tude lower than the ordinary ones, which significantly dimin- 
ishes the optically thin radiative force. Consequently, our mod- 
els showed that the contribution of less abundant hydrogen and 
helium isotopes 2 H and ^He to the radiative force is negligible. 

Similarly, we have also tested the contribution of lithium 
to the total radiative force. To this end, we have selected data 
on excitation and ionization of lithium ions from the Opacity 
Project database (Seaton [T987l Peach et al. 119881 Fernley et al. 
\TWj\ . The primordial lithium abundan ce N( \R)/N(lU) = 
4.15 x 10 -10 was taken from Coc et al. ©04). Numerical cal- 
culation showed that due to the extremely low lithium abun- 
dance its contribution to the radiative force can be also ne- 
glected. 

7.2. Expulsion of particular isotopes and 
abundance stratification 

Slightly different radiative force acting on different isotopes 
and different isotopic mass may cause radiative levitation of 
some isotopes and/or gravitation settling of others, the effect 
also referred to as the light induced drift (Aret & Sapar l200"2l . 
Thus, atmospheres of first generation stars may be chemically 
stratified. However, the discussion of these issues is beyond 
the scope of the present paper which is mainly intended to the 
study of an outflow from first generation stars. There appears 



a more interesting question in this context. Although the radi- 
ation force due to the less abundant isotopes is very small and 
has only marginal influence on the stellar atmospheres, it might 
be sufficient to expel these elements out from the stellar surface 
into the interstellar medium. 

However, numerical calculations showed that this is prob- 
ably not the case. Similar considerations already discussed in 
Sect.[6]for ordinary hydrogen and helium are valid also for deu- 
terium and the jHe isotope. Radiative force acting on the deu- 
terium is too low to allow a pure deuterium wind. Moreover, the 
radiative force due to free electrons is lower in ^He wind due 
to a lower number of free electrons per one nucleon compared 
to the H-He wind. 

Similarly, numerical calculations showed that pure lithium 
wind is also unlikely. Although the radiative force acting on on 
lithium may be higher than the gravity force in deeper layers 
of some models, in the outer regions is always too low to drive 
the lithium wind. 

8. Influence of line transitions on accretion 
physics 

Up to now, the initial mass function of first stars is not known 
(see Bromm & Larson 2004 1. One of the problems with the de- 
termination of initial mass function of first stars is connected 
with the process of termination of accretion onto a protostar. 
Dust grains are not present there, so the radiation pressure on 
dust grains, which operates in contemporary stars is not appar- 
ently available there. Thus, other mechanisms may be impor- 
tant for metal-free stars (e.g. Omukai & Inutsuka 2002 Tan & 
McKee l2004l 

Here we do not intend to study accretion physics in detail, 
but only to depict such processes which may be important for 
termination of accretion, namely we would like to test whether 
line transitions due to hydrogen and helium are able to influ- 
ence the accretion. 

Thus, we selected a fixed velocity law which can be ob- 
tained from the solution of the momentum equation for isother- 
mal (T = T w i n d) spherical accretion (Mihalas ll978l l 

/ a 2 \ dv r 2a 2 

\ Vr ! = 9- (33) 

\ v r J dr r 

The boundary condition of Eq. J33i was artificially selected to 
be v r (W0R*) = \J~\ a an d the density structure was calculated 
from the continuity equation 

4:Trpv r r 2 — M — const., (34) 

where M is the accretion rate. We selected three dif- 
ferent values of accretion rates, namely 1O -3 M0 year -1 , 
1O _5 M0 year -1 , 10 -7 M Q year -1 . The highest value of accre- 
tion rate roughly corresponds to the critical accretion rate over 
which protostellar core reaches the Eddington limit and no su- 
permassive star forms (Omukai & Palla 2003 1. 

For the velocity and density structure mentioned above we 
solved statistical equilibrium equations for hydrogen and he- 
lium inflow. We neglected any influence of possible shock at 
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the stellar surface and assumed that the protostar is in hydro- 
static equilibrium. In such a case for our numerical calcula- 
tions we may use the same stellar parameters as we used for 
wind tests, namely the stellar parameters taken from Kudritzki 
(2002 1. Our numerical calculations showed that the sum of the 
radiative force in the Sobolev approximation calculated after 
Eq. (I22t and the radiative force acting on free electrons can to- 
gether exceed the gravitational force only in stars close to the 
Eddington limit and only for relatively small accretion rates (of 
the order 10 _5 Mq year -1 and smaller). 

The principal difference between possible stellar wind and 
accretion flow is that important regions of accretion occur at 
much larger distances from the star than the acceleration re- 
gion of the stellar wind. Since the ionizing stellar flux at larger 
distances is lower, the inflow (especially helium) is less ionized 
and the line radiative force due to He II becomes important rel- 
ative to the gravity force. 

Clearly, the importance of the radiative force rises for more 
massive stars closer to the Eddington limit. Thus, there is a pos- 
sibility that for even more massive stars than those considered 
here the radiative force both due to line transitions and free 
electrons is able to overcome gravity and set the initial mass 
function of supermassive stars. These considerations have to 
be, however, tested against more advanced accretion models. 

9. Discussion of simplifying assumptions 

9. 1. Radiative transfer 

There are several problems which can influence our final re- 
sults. Probably the most serious one is the simplified radiative 
transfer in lines. This can influence the calculated occupation 
numbers and hence also the radiative force. However, this is 
probably not the case, since the results obtained are relatively 
robust, i.e. in many cases even the optically thin radiative force 
is not sufficient to drive a stellar wind. 

9.2. Level dissolution 

There is another more subtle problem which is related to the 
neglect of the so-called level dissolution (Hummer & Mihalas 
I1988> . It is the statistical treatment of the fact that the highest 
atomic levels cease to exist under the influence of surrounding 
particles. Usually this problem is treated approximately using 
the lowering of ionization potential. However, it can be treated 
in a consistent way using the so-called occupation probability 
formalism (Hubeny et al. U994l Kubat ll997b . Possibly, neglect 
of level dissolution can influence the calculated radiative force, 
specially in the case of hydrogen. However, test calculations 
with different number of hydrogen and helium lines showed 
this is not the case. Although total number densities of individ- 
ual ionization states may differ, number densities of lower ex- 
cited levels (which are important for line driving) are nearly the 
same. Since there is a relatively low radiative flux in the wave- 
length regions of lines originating from higher excited levels, 
their contribution to the radiative force is negligible (see also 
Sect. l2.2t . Thus, we conclude that the influence of level disso- 
lution on our results is only marginal. 



9.3. Pulsations 

Since we assume stationary flow here, the effect of possible 
stellar pulsations is neglected. According to numerical simu- 
lations of Baraffe et al. ( 12001 > . primordial stars may actually 
lose mass via pulsational instability. However, such effect is 
connected mostly with the stellar interior, so it is beyond the 
scope of the present paper. 



10. Conclusions 

10.1. Stars with chemical composition given by 
primordial nucleosynthesis 

We have studied the dynamical aspects of a possible circum- 
stellar envelope around massive first stars using our NLTE 
models. Although many of studied stars are very close to the 
Eddington limit and the maximum possible radiative force (cal- 
culated assuming that all lines are optically thin) is high in 
some cases, the stellar winds driven by line transitions are un- 
likely, because the line radiative force is very small. The influ- 
ence of other isotopes and elements, which had been produced 
during primordial nucleosynthesis, was found to be marginal. 
Also the effect of the continuum radiative force is negligible. 

We have also tested whether first stars can expel individual 
elements, either the dominant ones (i.e. hydrogen and helium) 
or the trace ones (i.e. hydrogen and helium isotopes or lithium). 
Calculations showed that the only theoretically possible chem- 
ically homogeneous outflow is that of hydrogen (for stars very 
close to the Eddington limit with T > 0.859), since other ele- 
ments yield a small number of free electrons per nucleon (and 
the line-radiative force is always negligible). However, even the 
case of a pure hydrogen outflow is unlikely since the radiative 
force in such a wind would be predominantly due to the free 
electrons. Hence, another process which would be able to sep- 
arate hydrogen and helium atoms in the atmosphere should ex- 
ists (e.g. gravitational settling) to enable a pure hydrogen wind. 
Without such a process a pure hydrogen wind is not possible. 

Finally, we have tested whether the radiative force due to 
line transitions is able to terminate the initial accretion onto first 
stars. We have found that this is only the case of stars relatively 
close to the Eddington limit and with relatively small accretion 
rates. Consequently, it is possible that also for even more mas- 
sive stars this process can influence the accretion physics and 
the initial mass-function of first stars. 

Our results do not imply that the first stars are completely 
without line driven winds during all of their life time. For ex- 
ample, cooler stars (which were not inlcuded into our present 
study) could have wind due to the hydrogen or helium lines. 
Moreover, if the accretion terminates, then as soon as sufficient 
amount of heavier elements is synthesised in the stellar core 
and transported to the surface layers, the line-driven wind is 
initiated (Kudritzki 2002). If this happens at a sufficiently low 
metallicity, such a stellar wind is multicomponent, allowing for 
element separation and frictional heating (Krticka et al. l2003l >. 
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10.2. General picture of stellar wind of massive 
low-metallicity stars 

From this study and previous ones (Kudritzki 2002 Krticka et 
al. 2003 1 a general picture of the properties of stellar wind of 
massive stars at low metallicities emerges. 

i Zero-metallicity stars with V < 0.859 likely do not have 
any stellar wind. If their atmospheres are relatively quiet 
(i.e. in the absence of violent pulsations or fast rotation) 
the chemical peculiarity in their atmospheres may de- 
velop. On the other hand, zero-metallicity H-Hc stars with 
r > 0.859 may have pure-hydrogen wind (probably very 
thin with mass-loss rates of order 10 -14 Mq year -1 , see 
Appendix iBl. driven only by the light scattering on free 
electrons. 

ii With increasing metallicity the line-driven wind may de- 
velop. For extremely low metallicities such a stellar wind 
is likely pure metallic (see Babel 1996 1. 

iii For higher metallicities also hydrogen and helium are ex- 
pelled from the atmosphere and multicomponent effects be- 
come important in such a stellar wind (Krticka et al. 2003 1. 
The multicomponent structure does not significantly influ- 
ence the wind mass-loss rate, which may be approximated 
by formulas obtained for one-component flow by Kudritzki 
(120021 . 

iv For higher metallicities, a line-driven wind which is similar 
to the stellar wind of present hot stars exists and its mass- 
loss rate was calculated by Kudritzki (2002i. 

Thus, in our opinion, the best choice for the evolutionary cal- 
culation now is to assume zero mass-loss rate for zero and low 
metallicity stars mentioned in items i-ii and to use Kudritzki 



( 2002 1 prescription for other stars. 
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Appendix A: Gravitational settling time scale 

Let us at least roughly estimate the characteristic time-scale td 
for gravitational settling of helium in the atmospheres of first 



stars. This is given by (Kippenhahn & Weigert 1990 p. 59) 



S 2 



(A.l) 



where S is the characteristic length for the density variation 
and the diffusion coefficient 



D 



kT 

377JH 



(A.2) 



The mean free path of helium in hydrogen plasma i we es- 
timate using the expression for Coulomb collision frequency 
(Burgers 1969 1. For typical parameters of the outer parts of the 
atmosphere (S = 10 11 cm, n « 10 8 cm" 3 , T w 30 000 K) we 
obtain mean-free path of £ ss 10 4 cm and corresponding char- 
acteristic diffusion time scale td ~ 10 4 years, which shows 



that the gravitational settling is really possible in the outer 
parts of atmospheres of first stars. On the other hand, in deeper 
atmospheric layers (where S — 10 11 cm, n ss 10 13 cm~ 3 , 
T w 50 000 K) we obtain mean-free path of £ « 1 cm and 
the diffusion time scale r D w 10 8 years, which highly exceeds 
the expected lifetime of massive first stars (Marigo et al. 2001 1, 
which is of order 10 6 years. Clearly, in the deeper layers the 
gravitational settling of helium probably does not take place. 

Anyway, the previous analysis enables us to roughly esti- 
mate the mass-loss rate due to the pure hydrogen wind, which 
is 

M H ~ AirRlSnm n T D l . (A.3) 

Note that since D ~ 1/n the estimated mass-loss rate does 
not depend on the atmospheric number density n, thus for 
both atmospheric points mentioned above we obtain Mh ~ 
1O _16 M0 year -1 . This value is close to a more exactly esti- 
mated upper limit of such hydrogen stellar wind, which is de- 
rived in AppendixlBl 

Appendix B: Pure hydrogen wind driven by free 
electrons 

The radiative acceleration due to the light scattering on free 
electrons may be higher than the gravitational acceleration 
in the pure hydrogen plasma for stars with T > 0.859 (see 
Eq. ( I32l i). Thus, if the hydrogen layer builds up at the surface 
of such star (e.g. due to the gravitational settling) then pure 
hydrogen wind may occur (cf. Malov [T974l for a description of 
stellar wind driven by Thomson scattering). In a stationary state 
such a stellar wind may be described (assuming spherical sym- 
metry) by hydrodynamical equations (Burgers 1969 see also 
Krticka & Kubat l200ll 

' {r 2 pHV r}1 ) = 0, (B.la) 



dr 



dv 



v rH - 



dr 



rad 
5e,H 



pu dr 



(4ph) 



-Kn a G(xHa), 
(B.lb) 



where indexes H and a stand for hydrogen and alpha particles, 
iVh an d Pu are the radial velocity and the density of completely 
ionized pure hydrogen plasma, = 2kT/mu (we assume 
that hydrogen and helium temperares are equal), the frictional 
parameter 



K Ha = n H n a k Ha = n H n a — In A, 



(B.2) 



and where the Chandrasekhar function G(x), defined in terms 
of the error function erf(x), is 



G(x) 



1 

2^2 



erf(x) 



2x 



: exp I 



and the argument of the Chandrasekhar function is 

v rH 



XHo 



(B.3) 



(B.4) 



where we have assumed that helium is static in the atmosphere 
and 

2 2kT (m H + m a ) 



a 



Ha 



mum a 



(B.5) 
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Inserting continuity equation Eq. JB.la> into Eq. IB. lbl and 
using the definition of Th (see Eq. J28» we obtain momentum 
equation in the form of 



we obtain equation for the Parker's type stellar wind (assum- 



dr 



GM 



do? 



(r H - 1) - ^ + - —k Ha G(x Ha ). 

dr r me 



(B.6) 



This equation has a critical point for u, H = an- To obtain a 
smooth solution at this critical point the critical condition 



GM 



(IW) 



dr 



24 



^k» a G ( *L 



(B.7) 



shall be fulfilled. Finally, the helium density in a static case is 
given by the equation of hydrostatic equilibrium in the form of 



5- (1 - r a ) — (a a p a ) H kn a G(x Ha ) = 0. 

r z p a dr m a 

(B.8) 

Eqs. (IB.6> - dB~8l (together with boundary condition which 
links hydrogen and helium density deep in the stellar atmo- 
sphere) can be used to calculate wind models of pure hydrogen 
flow. The mass-loss rate of such wind can be obtained by the 
requirement that the critical condition (IB.7I is fulfilled at the 
critical point. 

We do not aim to solve equations for pure hydrogen wind 
here, however we present only an estimate of upper limit of the 
mass-loss rate. To make its derivation more tractable, we ne- 
glect gas-pressure terms in Eqs. (IB.6> - dB~8l . The critical con- 
dition (1B.7I then yields the helium density at the critical point 



72— (!h- 1)t— 

r crit K H« 



G 



(B.9) 



Similarly, the gas pressure term in the helium hydrostatic equa- 
tion Eq. JB.8I > is always possitive, thus the upper limit for hy- 
drogen density at the critical point is 



it GM 



G 



a>ft a 



(B.10) 



With this we can derive upper limit for hydrogen mass-loss rate 
as 



M K <^ntv r rr 2 cA 



47ra H GM (1-T Q ) 



m H m, a 
kua 



G 



an 
ana 



(B.ll) 



Assuming q a — 2oh this can be rewritten as 



M H < 2xl(T 14 M Q year 



-if M 



T 



3/2 



(1 _ p ) 

V 100 m q y v K /" a 

(B.12) 

This shows that the mass-loss rate of hypothetical pure hydro- 
gen wind is very small. 

The wind velocity of such a wind outside the stellar atmo- 
sphere can be easily obtained from Eq. dB.6b assuming neg- 
ligible frictional force (G(xHa) -C 1), since in such a case 



ing isothermal wind, see Mihalas 119781 1. Note however, that 
the wind terminal velocity is influenced either by the radia- 
tive heating and cooling processes (which influence the sound 
speed an) or by hydrogen recombination, which decreases the 
radiative force on free electrons. 

Note however that the real situation will be likely more 
complex than that presented here. First, Eqs. JB.lt were de- 
rived assuming maxwellian distributions of velocities for all 
the components. While this may be true for hydrogen and he- 
lium, the velocity distribution of electrons may have two max- 
ima (one corresponding to electrons which move together with 
hydrogen and one corresponding to electrons which stay in the 
atmosphere). Possibly, also electric polarisation field, which 
was neglected in Eqs. JB. It . may play some role. Finally, the 
process described here may not be stationary. Anyway, it is not 
likely that pure hydrogen mass-loss rate is by orders of mag- 
nitude higher than that presented here, hence the influence of 
pure hydrogen mass-loss on the evolution of first stars is only 
marginal. 

For the calculation of the radiative force we have assumed 
that the stellar wind is optically thin for the Thomson scatter- 
ing, which seems to be a plausible approximation given the 
small mass-loss rate. 

It may seem surprising that we propose here that wind 
driven by scattering on free electrons may exist, whereas Malov 
i ll 9741 1 concluded that such wind may exist only in the case of 
a positive temperature gradient dT/dr > at the wind crit- 
ical point. His conclusion is based on critical point condition 
( IB.7I . which indeed cannot be fulfiled if dT/dr < and if 
we do not account for friction due to the helium atoms. Here, 
however, the friction due to the helium atoms enables to fulfill 
the critical point condition and acts in a similar way as positive 
temperature gradient in Malov's J1974I paper. 
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T40, dM/dt= 10" 6 M Q /year 
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Fig. 4. The plot of the relative difference (,g thin + gj ad - g) /g 
for model T40. Note that the points where this difference is zero 
match the points where the actual occupation number iVH e n,i 
is equal to N^\[ 1 relatively well. Thus, the simple conditions 
Eqs. J15I and Jl 81 for the minimum occupation number nec- 
essary to launch the stellar wind are able to predict relatively 
correctly the regions where the total radiative force is greater 
than gravity. 
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T50, dM/df=10~ b M„/year 



O) 
OS 
■o' 

+ 




T60, d/W/df = 10" b M„/year 



1 1.01 1.02 1.03 1.04 1.05 1.06 1.07 

r/R, 



T50, dimt= 10"' M„/year 



O) 

"3 

+ 




O) 
+ 

C 
x: 





-0.01 
-0.02 
-0.03 
-0.04 
-0.05 
-0.06 



1.01 1.02 1.03 1.04 1.05 1.06 1.07 

r/R, 

T50, d/W/df=10" 8 M„/year 



"C 

+ 




1 1.01 1.02 1.03 1.04 1.05 1.06 1.07 

r/R, 

T60,d/W/df= 10" 7 M„/year 





0.1 r- 




0.08 - 




0.06 




■ 




0.04 l 








0.02 - 


+ 




JZ 


L 








-0.02 : 




-0.04 : 




-0.06 - 
1 




r 




-0.01 




-0.02 : 




-0.03 




-0.04 : 


+ 


a 




'sz 


-0.05 - 




-0.06 b 




1 1.01 1.02 1.03 1.04 1.05 1.06 1.07 

r/R, 



-0.07 



r/R, 

T60, d/W/df= 10" 8 M Q /year 



1 1.01 1.02 1.03 1.04 1.05 1.06 1.07 

r/R, 



Fig. 5. The same as Fig.|4] however for model T50. 



Fig. 6. The same as Fig.0] however for model T60. 



